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The contribution of pseudoscalar pole diagrams to the hadronic light-by-light corrections of order 
a 3 Ep to the ground state hyperfine splitting in hydrogenic atom is calculated in the pseudoscalar 
pole terms approximation. The vector dominance model for the form factor of the transition of a 
pseudoscalar meson into two photons is used. With the account of new experimental data on the 
cross sections a(e + e~ — > p,uj — > 7r°(7;)7) in the SND and CMD-2 experiments some other hadronic 
corrections to the muonium hyperfine splitting are calculated. 
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The study of strong interaction contributions to the energy spectra of hydrogen-like systems and lepton anomalous 
magnetic moments (AMM) is considered lately as an important problem which has large practical meaning for testing 
the Standard Model Jl|,||. One of such contributions is determined by hadronic vacuum polarization (HVP) |,|. 
Presently there are at least two tasks where the HVP effects can be tested on the experiment. The measurement of 
muon anomalous magnetic moment in the new E821 experiment at the Brookhaven National Laboratory (BNL) was 
done with extremely high accuracy : 

a p = 116 592 023 (151) x KT U , (1) 

which allows to perform a check of the hadronic contribution in the Standard Model not only at one-loop level, but 
also with two-loop accuracy. In a Los Alamos experiment aimed at measuring the hyperfine splitting (HFS) of the 
muonium ground state, the accuracy of the measurement reached a few tens of Hertz 

A^ HFS (Mu) = 4 463 302 765 (53)Hz. (2) 

This requires taking into account both higher order contributions in a (a is the fine structure constant) and the 
contribution from the HVP of different order |?]|32|,^|,[lO| to the muonium hyperfine splitting. Consistent calculation of 
the HVP contribution to the photon Green's function can be carried out on high level of the accuracy, because hadronic 
spectral function p h is connected with the cross section of e + e" annihilation into hadrons measured experimentally. 
A completely different situation occurs in the case of the hadronic light-by-light scattering contribution to the HFS 
in hydrogenic atoms or to the muon AMM [pl}]l4| (see diagrams in fig. 1,2). The contribution of such diagrams to 
the energy spectrum or lepton AMM is determined by the tensor of virtual 7*7* scattering containing eight structure 
functions. Reliable experimental data about these structure functions are not available at present. In the low energy 
region 7*7* interaction may be studied by means of effective field theory of interacting photons and hadrons f p"5| . 
The calculation of such hadronic contributions to muon AMM in the pion pole terms approximation was carried out 
in Ref. [^6) . As a result the difference between the theoretical value for the muon AMM and experimental data was 
reduced to about 1 standard deviation P,|l7|^9|. In this work we calculated the contribution of pseudoscalar pole 
terms to the muonium HFS Q. Corresponding diagrams, shown in fig. 1,2 describe a part of O(a) corrections to the 
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leading hadronic contribution in the muonium HFS. The investigation of higher order effects on a to the muonium 
ground state hyperfine splitting allows to determine the muon-to-electron mass ratio from experimentally measured 
value (2) 

— ^ = 206.768 279 8(43), 6 = 2 • 1(T 8 . (3) 

m c 

Such indirect determination of the muon electron mass ratio is almost six times more accurate than the best experi- 
mental value for m^/m c Further improvement of the accuracy in (3) is connected with the calculation of a new 
theoretical contributions to Ai^ s (Mu) of order 10 Hz, containing small parameters a 4 Ep ~ 13 Hz, a 3 j^Ep ~ 8 Hz, 

a 2 (m c /m At ) 2 Ep k 6 Hz which can be enhanced by the presence of logarithmic factors lna -1 = 4.92, lnfm^/m,,) = 
5.33. 




+ 




FIG. 1. Pseudoscalar pole terms giving rise the contribution of order a 3 ( m °™^ J Ef to the muonium HFS. 
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FIG. 2. Pseudoscalar pole terms giving rise the contribution of order a 3 ( mi ^ v ^ Ef to the muonium HFS. 



The effective vertex for the interaction of a ir° meson ( or other pseudoscalar mesons 77, r/) and virtual photons can 



be expressed by means of the transition form factor F„.o 7 * 7 » (k 2 , k|) 



V^(k!,k 2 ) 
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F 7r o 7 . 7 .(k 2 ,k 2 ) 
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where F^ = 0.0924 Gev is the pion decay constant, ki,k2 are the photon four momenta. Different model represen- 
tations for the transition form factor F 7r o 7 » 7 * (k 2 , k 2 ,) exist. All of them can be used for the numerical estimation of 
the corrections to the muonium HFS. Firstly consider corrections which are determined by the diagrams in fig.l. In 
this case the photons interacting with one particle (electron or muon) possess the momenta ki^ and (— ki.2). The 
Feynman amplitudes containing the photon momenta ki.2 and (— ki,2 + qi, 2 — Pi. 2) (pi,2 and qi.2 are four momenta 
of particles e _ , /i + in the initial and final states respectively) in one vertex (4) can be omitted since they determine 
the contribution to HFS of higher order in a. 

To construct the HFS part of the quasipotential in the e" 
states with the spin S=0 and S=l respectively: 



fi + system we used the operators projecting onto the 
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(5) 
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where e M is the muonium polarization vector for the state 3 Si, e = £^7^. To obtain the contribution of the Feynman 



amplitude M to the energy spectrum we must calculate the common trace: Tr 
photons in fig. 1 (a) give the following contribution to the muonium HFS: 



ft+Mft 



The diagrams with direct 
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where the Fermi energy Ep = 8/i (Za) /m e m M , /i = m c m Al /(m + m^) is reduced mass, n is the principal quantum 
number, 



fl(ki,k 2 ) = -2k° 2 k^(k!k 2 ) + k^^k 2 + (kxk 2 ) 2 ] + (k!k 2 )[kfk^ - (kik 2 ) 2 ]. 



(7) 



We make dimensionless the integration variables ki , k 2 using the mass of the vector meson My ■ Our function f i (ki , k 2 ) 
was derived after the calculation of the trace, summing over the Lorentz indices and averaging over the directions of 
the three momentum k for the states with total spin S=0 and S=l respectively: 



e Mra/3 k! Q k 2 «e A ^k! p k 2 w Tr 



7 A ( Pl - ki + m e )^H^(-p 2 + t 2 + m M ) 7 CT n^ 



(8) 



We employ the system FORM for this aim p3| . The contribution of the Feynman diagram in fig.l (b) to HFS with 
crossed photons has the form: 
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(9) 



Corrections of lower order in the ratio of particle masses to the muonium HFS are determined by the Feynman 
diagrams in fig. 2. They contain the factor of the same power in a, because the photon propagator 1/k 2 enters the 
amplitudes together with two factors k in the numerator. One degree of k is related to the anomalous part of the 
electron transition current u(qi)^ 4 -cr A1,y k ly u(pi) and the other enters one of the vertices of the 7r°7*7* interaction. 
This does necessarily mean that the muon transition current expressed as 
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u(q 2 )r^u(p 2 ) = u(q 2 ) 
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(10) 



has the contributions from diagrams in fig. 2 only to the magnetic form factor F m , while F = 0. The sum of 
contributions coming from the diagrams in fig. 2 (a) and (b) can be written in the unified form if we use the symmetry 



properties of integral function under changing the variables ki 



-k 2 : 
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The contribution of the diagram in fig. 2 (c) which has different structure of the pion propagator and the transition 
form factors F 7r o 7 » 7 * (k 2 , k 2 ) can be presented as follows: 
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where 
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Taking into account that the energy corrections (11) and (13) have the structure Epa M we conclude that hadronic 
contributions to muon AMM from expressions (11) and (13) coincide with the relating results of Ref. To extract 
the anomalous part of muon electromagnetic current containing the form factor F 2 , the authors of Ref. [ fl6| used the 
projection operator of special type. The form factor F^ 7 « 7 »(k 2 ,k|) describing the conversion of two virtual photons 
into 7r meson should be calculated on the basis of strong interaction theory. The asymptotic behaviour of this form 
factor for large photon virtualities k 2 2 was studied using the methods of perturbative and nonperturbative quantum 
chromodynamics ]24| , |25| ] . But the calculation of the corrections to the muonium HFS requires the proper description of 
the form factor F 7r o 7 » 7 » (k 2 , k 2 ) at small and intermediate photon virtualities k 2 , k 2 around 1 Gev 2 , which are dominant 
in the integrals entering AEp FS . To evaluate the obtained contributions to the muonium HFS it is convenient to use 
the representation of the form factor F„.o 7 . 7 , (k 2 , k 2 ) based on the vector dominance model (VDM): 



(k?,k 2 ) = 



1 



1 
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(15) 



where in the case of the tt contribution Ay = My = M p =0.7693 Gev is the mass of p meson. This double VDM 
ansatz which was used as a theoretical input in the CLEO data analysis J2(|, behaves as l/(k 2 k 2 ) for asympotically 
large k 2 , k|, contradicting predictions of pQCD. Nevertheless for reasonably small momenta k 2 , k 2 , this model agrees 
with the experimental data p6| , so we used it to obtain the numerical evaluation of relevant hadronic contributions. 
We performed numerical integration in the expressions AEp FS using (15). As a result we obtained the summary 
contribution to the muonium HFS (AE^ FS + AE 2 iFS )(7r°) coming from the diagrams in fig.l which is equal 0.001 Hz. 
Respective result for the sum of the diagrams in fig.2 (AE" FS + AEJ fs )(tt°) = 2.515 Hz. The employment of the 
VDM for the form factor F 7r o 7 « 7 » (k 2 , k 2 ) makes possible analytical calculation of the integrals (11), (13), which was 
done in Ref. p8| by means of the expansion of integral functions in two parameters S = (m 2 — m 2 )/m 2 and m 2 /My. 
The particle interaction amplitudes shown in fig.l have additional small factor m e /M p in the comparison with the 
amplitudes in fig.2. So, corresponding numerical values are very different. We note also that the numerical value of 
the contribution from the diagrams in fig.2 to the muonium HFS a c Ep which is determined by the electron anomalous 
magnetic moment is extremely small: 10 -4 Hz. 
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FIG. 3. Fit of the SND cross section cr(e+e~ -> tt°7) 

The contributions of other pseudoscalar mesons r\ and rf were estimated in a similar way. To describe couplings of 
r\ and rf with virtual photons we used the form factor representation (4) where the pseudoscalar decay constants F^ 
and F v i are related to the two photon partial width of the resonance [^6| : 

64tt 3 r(i7, rf -> 77) ' V ; 

The CLEO data lead to the following values of F P ||]: F„ = 0.0975 Gev, F v > = 0.0744 Gev. Their fit for transition 
form factors with a function defined by eq.(15) gives the following results for the parameter Ay: A, ; = 0.774 Gev, 
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A V ' = 0.859 Gev. The results of numerical integrations for these two pseudoscalar mesons and total contributions to 
the muonium HFS are presented in the table. Note that from the theoretical viewpoint F,,, entering in Eq. (16) 
should be considered as effective decay constants due to the r) — r]' mixing [^7j . We considered here the pseudoscalar 
pole terms contributions of order a 3 Ep to the muonium HFS presented in Fig. 1, 2. Numerical values of corresponding 
contributions are very different. The contributions of diagrams in Fig. 2 can be taken through the muon anomalous 
magnetic moment since the value of the muonium HFS is proportional to the product of the electron and muon total 
magnetic moments. So the new contributions to the muonium HFS come from the diagrams in Fig. 1 and are given 
in the first column of the table. 




FIG. 4. Fit of the CMD-2 cross section cr(e + e" -» 777) 




TABLE I. Contributions of pseudoscalar mesons to light-by-light hadronic interaction in the muonium HFS. 



p 


(AE HFS + AE? FS ), Hz 


(AEf FS + AE HFS ), Hz 




1.1 x 10~ 3 


2.515 


V 


0.2 x 10" J 


0.586 


rf 


0.1 x 10-" 


0.552 


Total sum 


1.4 x 10~ a 


3.653 
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The contributions of the annihilation processes e + e~ — > p,u> — » ir°-f,r)j to the HVP corrections to the muonium 
HFS have the same order of the magnitude at energies y/s < 1 Gev. Recently new measurements of the cross sections 
of these reactions were carried out at SND and CMD-2 (Novosibirsk). The obtained experimental data for the cross 
sections a(e + e~ — > tt°j) and <r(e + e _ — > 777) and their theoretical approximations are shown in fig. 3, 4. Reactions of 
this type were taken into account effectively in our previous calculations 0j of the HVP corrections to the muonium 
HFS by using the Breit-Wigner parameterization for the resonance to. There are three important decay channels 
of the resonance u>, which contribute to the total decay rate with the per cent accuracy: u> — > ir + ir~ir° (88.8%), 
uj — > 7r°7 (8.5%), ojtt + tt~ (2.21%) pq |. The appearance of new experimental data for the a(e + e~ — > u — > 7r + 7r~7r ), 
cr(e + e _ — ► p, u> — ► 7r°(?7)7) [p9]-|3"lf allows to separate hadronic contributions of the processes to the muonium HFS and 
to perform the calculation with higher accuracy. Theoretical cross sections for these processes in the region of p, u> 
mesons were presented in the form p9|--pl|: 



f (s) = Fn °^ |A p07r o 7 (s) + A W7r o 7 (s) + A 07r o 7 (s)| 2 , (17) 



^ 7 (s) = ^|V4^C,+ Yl A v.7l 2 > (18) 

V— p,LU,4>,p' 



^(^^^^•lA^+e-A^ + Abgl 2 , (19) 



Av *W B > " m 2 ,-s-iVirv(s)V CTV " 0W7 F7m v T' ( } 

where F phase space factors of the final states, ervu-o(?j)7 is the cross section of the process 

e + e~ — » V — » 7r°(?7)7 at ^/s = my, my is the resonance mass, Ty(s) is its decay width at the squared cm. en- 
ergy s, Ty = ry(m v ). Three pion decay amplitudes A w , A^, have the parameterization similar to (20) F 37r (s) is 
a smooth function which describes the dynamics of V — > p7r — > 7r + 7r _ 7r° decay |3^]. The form of the functions entering 
the expressions (17)-(20) and the values of numerous parameters of the cross sections (17) - (19) were taken also from 
Ref. p9| The theoretical fit (19) of the experimental data is shown in fig. 5. To obtain the contributions to the 
muonium HFS we performed numerical integration of the cross sections (17) - (19) as in the previous work [Q in the 
energy region 0.6 < t/s < 0.81. As a result we have found the following values of the corrections: AE^ V p(p, u> — * 7r°7) 
= 1.33 ±0.03 Hz, AEg v |(p,w -> 777) = 0.042 ± 0.003 Hz, AEg v |(w -> 3?r) = 10.82 ± 0.33 Hz. So we obtain that our 
previous result for the ui contribution from 12.45 Hz must be changed on 12.19 Hz. 

In conclusion we can point out that among other hadronic contributions of order a 3 Ep there are one-loop corrections 
which contain in parallel with hadronic loop also the electronic or muonic vacuum polarization (EVP and MVP 
respectively). The expressions for these contributions can be obtained using the skeleton integrals for two photon 
exchange diagrams ]33| , p0] | and standard modification of the photon propagator due to vacuum polarization and 
hadronic vacuum polarization. Without going into further details we present for these corrections the numerical 
values: AE^f HVP - 3.7 Hz, AEg™ (HVP = L1 Hz . 

We have calculated all hadronic light- by- light corrections of order a 3 Ep in the pseudoscalar pole terms approxima- 
tion and main one-loop contributions containing electron (muon) and hadronic vacuum polarization of the same order 
in a to the muonium HFS. Combining all results presented here we find that the theoretical value of the hadronic 
contribution to the muonium ground state HFS should be increased by 8.3 Hz. 
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